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Abstract. We study the influence of the assumption behind the 
use of the inverse Tully-Fisher relation: that there should be no 
observational cutoffs in the TF parameter log Vm- It is noted 
how lower and upper cutoffs would be seen in a \ogH vs. 
"normalized distance" diagram. Analytical expressions, under 
the simplifying assumption of a normal distribution and the use 
of the correct TF slope, are derived for the resulting biases, 
especially the average bias which log Vm cutoffs produce in the 
derived value of Hq. This bias is shown to be relatively weak, 
and as such cannot explain the large differences in the reported 
values of Hq derived from direct and inverse TF relations. 

Some problems of slope and calibration are shown to be 
more serious. In particular, one consequence of fitting through 
the calibrators either the slope relevant for field galaxies or the 
steeper slope followed by calibrators is that the derived value 
of the Hubble constant comes to depend on the nature of the 
calibrator sample. If the calibrator sample is not representative 
of the cosmic distribution of log Vm, large errors in the derived 
value of Hq are possible. Analytical expressions are given for 
this error that we term the calibrator sample bias. 

Key words: galaxies: spiral - galaxies: distances and redshifts 
- cosmology: distance scale 



1. Introduction 

There have been sincere hopes that the inverse Tully-Fisher 
(TF) relation could overcome the distance dependent selection 
bias, the so-called Malmquist bias of the 2nd kind (Teerikorpi 



1997), influencing the determination of the value of the Hubble 



constant. First, Schechter (198C) pointed out that the inverse 
relation of the form 



p = a'M + b' 



(1) 



where M is absolute magnitude (or linear diameter) and p is an 
observable parameter not affected by selection effects (e.g. the 
maximum rotation velocity log Vm in the TF relation), can be 
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derived even from magnitude limited samples in an unbiased 
manner. Secondly, Teerikorpi (1984, or T84) showed explic- 
itly that the inverse relation gives unbiased average distances 
for a sample of galaxies (if p is not affected by selection), 



a result later confirmed by Tully (1988) via simulations (see 



also Hendry & Simmons 1994). Thirdly, Ekholm & Teerikorpi 
(1997, or ET97| ) showed that with the inverse relation, one does 
not necessarily need a volume-limited calibrator sample, which 
on the contrary is quite critical for the use of the direct relation: 



M = ap + b 



(2) 



The direct TF relation has been succesfully applied for deriva- 
tion of Hq by the method of normalized distances that was 
first used by Bottinelli et al. ( |1986| ). Theureau et al. ( |l997b| ) 
have used an improved version of this method in an analysis 
of a sample of 5171 spiral galaxies that leads to the value of 
Hq ss 55 ± 5 kms -1 Mpc -1 . In this method, which is rather 



similar to the Spaenhauer diagram method of Sandage (1994), 
one identifies the so-called unbiased plateau in the (log H) vs. 
normalized distance diagram for a magnitude or angular diam- 
eter limited sample, from which H is determined without the 
need to make further model dependent correction due to se- 
lection. Normally, the unbiased plateau contains only 10-20 
percent of the total sample. That one can use the whole sample 
without bias correction is a remarkable advantage of the inverse 
relation approach. 

However, several problems balance the advantages of the 
inverse relation (Fouque et al. 1990, Teerikorpi 1990 Sandage 
et al. 1995, ET97), and one has not yet been able to use it re- 



liably for an independent determination of Hq (see Theureau 
1998| ). In the present paper, we study two 
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important points which must be clarified before a successful 
use of the inverse TF relation: 

1 . How would possible cutoffs in log Vm be seen in the data? 

2. How large an influence do upper and lower cutoffs in 
logVM have on derived average distances (and hence on 

Hq)1 
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3. How does the nature of the calibrator sample influence the 
derived average distances if the inverse TF slopes for cali- 
brators and field galaxies differ? 

As regards the first point, we emphasize that our approach 
follows the lines which we have adopted and found useful with 
the direct TF relation: to detect and overcome the bias without 
detailed modeling, by studying how the cutoffs would appear 
in kinematical applications (c.f. the increase of (H) as a func- 
tion of true (or normalized) distance, when the direct TF rela- 
tion is used). Another, mor e tech nical approach to the direct TF 
relation is that of Willick (1994), based on corrections to indi- 



vidual galaxies within an iterative scheme. In principle, a good 
knowledge of the log Vm selection function could also permit 
Willick' s method to work for the inverse TF relation, in order 
to derive the slope as it would be without the selection. 

Our approach necessarily cuts a way a part of the sample. 
However, as shown by Ekholm et al. ( 1998 ) for the KLUN sam- 
ple, the remaining subsample suitable for the determination of 
Ho with the inverse TF method is still much larger than the 
unbiased subsample for the direct TF relation. 

The observational motivation for studying question 2) 
comes from Theureau et al. (1998), who discuss the detec- 



tion rate of 21cm line profiles at the Nancay radio telescope, 
in connection with the angular size limited KLUN galaxy sam- 
ple. The detection rate is relatively high, 86 percent for galax- 
ies with a previous redshift measurement and 61 percent for 
those with unknown z. However, more important than the av- 
erage detection rate, is how detection depends on p, i.e. the 
selection function S(p). Unfo rtunat ely this is difficult to derive 
from the raw data. Theureau ( 1997 ) notes that one expects loss 
of galaxies with either a narrow or broad 21cm line width. Nar- 
row profiles are difficult to detect among the noise, while broad 
profiles, being low, also tend to be missed. 

Hence, as a first approximation we suppose that the distri- 
bution of log Vm is affected by sharp lower and upper cutoffs, 
p\ and p u , although the selection S(p) more probably depends 
on various galactic parameters in a complicated manner. How 
such cutoffs influence, e.g. the value of Hq as determined by 
the inverse relation method, has not been previously quantita- 
tively discussed. 

We study the above questions using normal distributions, 
which allow analytic expressions. This forms a natural sequel 
to our previous discussions of the TF relations, both direct and 
inverse, where the assumption of Gaussianity has been adopted 
(Teerikorpi 1984, |1987|, |1990[ |1993[). It is worth noting that 



though we occasionally refer to the KLUN sample (see e.g. 
Theureau et al. 1997a), the theoretical ideas concerning the 



above mentioned points 1-3 have general validity. As an ex- 



tensive sequel to the present work, Ekholm et al. (1998) have 
investigated whether the large value of Hq, derived with a 
straightforward application of the inverse TF relation for the 
KLUN sample, could be due to the errors arising from cutoff 
and/or calibrator sample biases. They show that the calibrator 
sample bias is more important. 
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Fig. 1. p vs. log of apparent diameter, log D25 diagram for dis- 
tant cluster and calibrator galaxies. Upper cutoff at p (dotted 
line) makes the inverse TF distance to the cluster too small 
when we force the correct slope, followed by both the cali- 
brators and the cluster (dashed lines), through the cluster. The 
solid vertical line shows the angular diameter limit, the dashed 
vertical line is the 'unbiased plateau' limit for the distant clus- 
ter. Cluster members that are not in the observed sample are 
marked with small symbols. 

In this paper we also aim at a better general understanding 
of the inverse relation as a distance indicator, and present some 
arguments in a heuristic manner, since the method has been 
justly criticised as difficult to visualize, in comparison with the 
"more natural" direct relation approach. 

This paper is structured as follows: In Sect. 2 we show qual- 
itatively how cutoffs in log Vm influence the distance determi- 
nation. In Sect. 3 the method of normalized distance d n is in- 
troduced for the inverse TF relation. In Sect. 4, the behaviour 
in the logiJ vs. d n diagram is calculated analytically. In Sect. 
5, we derive an analytical expression for the average bias in 
logiJ when there is an upper or lower cutoff in log Vm, and 
give some examples in Sect. 6. The fundamental problems of 
slope and calibration are discussed in Sects. 7 and 8. Section 9 
contains concluding remarks, emphasizing the essential points 
and consequences of the present study. 

A note on the nomenclature: We use capital D's for the di- 
ameter of a galaxy, D25 is the usual apparent 25th-mag isopho- 
tal diameter ([D25] =0.1 arcmin), D is the linear diameter 
([D] = kpc). In Sect. 4, to make the formulae more readable, 
we use T> for the log D's. Small d's are used for inferred dis- 
tances, d n is the normalized and dkin the kinematical distance. 
It is assumed that kinematical distances have been calculated 
from a realistic velocity field model, giving reliable relative dis- 
tances (in practice from a Virgo-centric infall model as e.g. in 
Theureau et al. 1997b). The true distance is denoted by r. 



2. Influence of cutoffs in log Vm: schematic treatment 

Let us assume that the selection may be described as cutoffs on 
the wings of the symmetric p distribution. To illustrate, let us 
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consider the simple case when there is only an upper cutoff p u , 
so that all p > p u are missed, while all p < p u are measured. 
The situation is sketched in Fig. 1. The calibrator sample (at a 
fixed distance) is on the right, while a distant cluster of galax- 
ies is on the left. The p u -limit is indicated by a dotted line. 
The correct regression lines, with the same slopes, are shown. 
Note that unaware of the p selection, one would derive from 
the cluster galaxies (or more generally from field galaxies) a 
too shallow slope. For the moment we simply assume that the 
correct slope is known. We return to the problem of slope in 
Sect. 7, and this issue will be discussed elsewhere (Ekholm et 



al. |1998j ) in connection with real data. 

In inverse TF method (iTF) for deriving distances, one ba- 
sically shifts one of the regression lines over the other, and the 
required shift of log D25 gives the (logarithmic) ratio of cluster 
to calibrator distance. 

In practice, one calculates for each cluster galaxy the iTF 
"distance" d c \ using the calibrating relation p = a' log D + b', 
so that log (<i c i /4ai) = log D(p) - \ogD 25 - Averaging over 
all cluster galaxies gives the unbiased distance. This is equiv- 
alent to "shifting the regression lines". The iTF "distance" has 
a systematic error depending on p, t houg h one does not know 
how large the error is (c.f. Sect. 5 in T84 ). 

In the depicted situation, even if one knows the correct 
slope, the attempt to force-fit it on the cluster data will shift 
the line towards larger logics, because the upper cutoff has 
eliminated more galaxies from above the "correct" (unbiased) 
regression line than from below it. Hence, the inevitable result 
of an upper cutoff is that the iTF distance to the cluster will be 
underestimated. 

Now, because of the limiting angular size, not all galax- 
ies of the cluster are included. So one should add to the di- 
agram a vertical line, only galaxies to the right of which are 
large enough to enter the sample. Note that this cutoff does not 
bias the distance indicator if there is no p-cutoff. 

If one moves the cluster closer (to the right), the angular 
size limit cuts away less of the cluster galaxies. Then the influ- 
ence of the p u -cutoff on the derived distances becomes weaker. 
In terms of the calculated Hubble constant (actually up to an 
unknown factor, because b' /a' is still unknown), one should 
see in a \ogH vs. dki n diagram first at small kinematical dis- 
tances (proportional to the true ones) a constant "plateau" when 
the angular size limit keeps out of the cluster. Then there is an 
increase due to the larger influence of the p u -cutoff. 

3. "Unbiased plateau" for the inverse relation 

The plateau in the logH vs. dki n diagram is not an unbiased 
plateau; it still gives too large a value of (logH), because of 
a residual effect of the p u -limit. The bias can be reduced by 
applying an upper cutoff to the cluster sample in logI?25, as 
illustrated by the dashed vertical line in Fig. 1. The galaxies 
between the upper and lower log D25 cutoffs now provide the 
best iTF distance estimate. 

The position of the critical righthand side vertical line in 
Fig. 1, which cuts away the "bad" logZ?25 interval, depends 



on the distance of the cluster (compare the cluster and the cal- 
ibrators). However, its position is practically constant on the 
linear size scale, independently of the cluster distance. This is 
clearly seen if one shifts in Fig. 1 the cluster galaxies to the 
right, i.e. decreases their distance until it coincides with the 
distance of the calibrator cluster. Because the TF relations are 
essentially the same, the critical vertical lines will almost coin- 
cide. Hence, in the general case of a field galaxy sample, one 
may shift all galaxies independent of their distances, to one di- 
agram of log Vm vs - logI?25 + logdkin, where the abscissa is 
up to an unknown constant (as far as Hq is unknown) the same 
as linear diameter, log D. If one now accepts galaxies left of 
the (common) critical line only, an unbiased average distance 
is derived. 

One might term log Z? 2 5 + log d k i n (or d n = D25 ■ d kin ) 
the normalized distance, analogously to th e case of the direct 
TF normalized distance introduced in 184 Then in the log H 
vs. log d n diagram, using the correct slope a', there is at small 
normalized distances an unbiased plateau, after which (log H) 
starts to grow. 

The resulting unbiased plateau should be horizontal, if the 



slope a' is correct. ET97 pointed out a simple device of decid- 
ing which is the correct inverse slope: in the log H vs. dkin dia- 
gram the run of points should be horizontal. This was, however, 
made under the premise of no p-selection. Now we understand 
that instead of dkin one should use d n = £>25dkin on the ab- 
scissa, in order to identify the cutoff d n . Only galaxies below 
this cutoff are expected to show a horizontal run on the log H 
vs. dkin diagram, if a' is correct. Figure 2 shows how in the 
more general case, the lower cutoff caus es too small (log H) 
below the unbiased plateau. Ekholm et al. 1998 will discuss the 
evidence for upper and lower cutoffs in log Vm in the KLUN 
sample and how the described method can in practice be im- 
plemented. 

Though in this manner one may detect sharp cutoffs in 
log VM and derive an unbiased value of Hq without detailed 
modelling, it is important to know how large an influence the 
cutoffs would have on Hq, if their presence were ignored. We 
discuss this in Sects. 4-6. 



4. Analytic calculation of (log H) vs. d n 

Assuming Gaussian distributions and sharp cutoffs atpi andp u , 
it is easy to calculate the expected run of (logH) vs. normal- 
ized distance d n . For simplicity, we write the inverse relation 
as 



p = a'V + b' 



(3) 



We use T> = log D\ m as shorthand notation for the logarithm 
of the linear diameter. At any fixed T> the dispersion of p is a v . 

Now the task is to calculate the bias in log H at eac h T>. 
This is done starting from the same initial formula (28) of T84 
which was used to show the absence of bias. However, in this 
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Fig. 2. Effect of the p-cutoffs on the value of (log H) at differ- 



case one must add the selection function S(p): 

/~„(X>-2>(p))S(p)3(p)dp 



(bias)p = 



const • exp [— (p 



(4) 



where $(p) = const • exp [— (p - p(D)) 2 / (2cr|)] and V(p) = 
[p-b')/a'. 

This formula gives for the sample galaxies with true (log) 
diameter T>, the (log) error in the diameter deduced by the in- 
verse relation (using T)(p)). It is directly reflected to the de- 
rived value of (\ogH) at constant T> (i.e., at fixed normalized 
distance d n ). 

With the present choice of the selection function (sharp up- 
per and lower cutoffs), and after a change of variable, 



p-a'V-b 1 

Eq. (4) becomes 

(bias)p = -yj2/-K(o p /a')B(x\,Xu) 



(5) 



(6) 



Here the term B(x\, x u ) can be evaluated in terms of the error 
function erf (x) = (2/y / 7r) J Q exp (~t 2 )dt: 



B(x h x u ) 



exp (— xf ) — exp (— x^) 
erf(x u ) — erf(xi) 



(7) 



where x\ and x u correspond to the cutoffs p\ and p u via Eq. (5). 

Inspection of the formula confirms the qualitative conclu- 
sions of Sect. 2. Putting x\ = — oo, so that only the upper cutoff 
remains (and recalling that erf (— oo) = —1), the bias is seen to 
be positive: the diameters are underestimated, leading to under- 
estimated distances, hence to overestimated \ogH. Note that 
this calculation does not tell us the average bias in the sample. 

5. Calculation of the average bias 

In order to have an idea of how large an influence the p-cutoffs 
have on the value of (log H) calculated from the whole sample 




Fig. 3. A volume limited subsample of a synthetic data-set is 
shown as crosses, and the diameter limited subsample of the 
same parent data-set is marked with circles. Direct TF relations 
of these subsamples, £>dir(p) and 2?Ma,dir(p), are parallel and 
separated vertically by (3 — a) ■ In 10 • cr|, (Eq. 11). The in- 
verse relation T>- lm {p) (Eq. 12), that under ideal conditions is 
the same for both subsamples, intersects the direct regression 
lines at V = V and V = V + (3 - a) 



In 10- o-f,. 



using a (correct) inverse relation, we derive the average bias for 
a special case of radial space density behaviour: 



n(r) = kr 



(8) 



where r is the radial distance from our Galaxy and a = cor- 
responds to homogeneity. It is easy to understand that the aver- 
age bias depends on how the observed galaxies are distributed 
along the d n -axis, which must depend on how galaxies are dis- 
tributed in space. 

Rather than trying to start with averaging the derived 
(bias)d n over d n (Eq. 6), we use another route, based on the 
fact that the bivariate distribution F(d n ,p) is known for an an- 
gular diameter limited sample, when space density behaves as 
Eq. (8). We assume for simplicity that the cosmic distribution 
of p (as well as of T>) is normal, which is a rough approxima- 
tion for fixed Hubble types, see e.g. Fig. 1 in Theureau et al. 



(1997a) 



As illustrated by Fig. 1 in T84 for the case of magnitude 
TF relation, in magnitude and volume limited samples the di- 
rect regression lines are parallel and separated by 1.382er^ f 
(when a = 0). The separation happens so that the volume- 
limited line (carrying with it the bivariate distribution) glides 
along the unchanged inverse line which goes through the av- 
erage M of both the volume ( = Mo) and magnitude limited 
( = Mq — 1.382(7^) bivariate distributions. Figure 3 gives a 
similar diagram for the diameter TF relation, showing the three 
kinds of regression lines for a synthetic sample of field galax- 
ies. 

In the present calculation we need to know how the diam- 
eter direct r elatio n is constr ucted f rom the parameters defining 
the iTF (c.f. T84 Teerikorpi 1993, where this was given for the 
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magnitude TF relation): 

X>dir(p) = (a'/a 2 p )a 2 Vp ■ p + [T> /o* - a'b'/a^a^ (9) 
where 

V4 P = 1/4 + VK/a') 2 - (10) 
The Malmquist shifted (angular size limited) direct relation is 

£>Ma,dir(p) = Pdir(p) + (3 - a) • In 10 • a% v (11) 

and the inverse relation is 

-D im (p) = {l/a')p-b'/a'. (12) 

One can express the average bias at a fixed p using the in- 
verse and (Malmquist-shifted) direct relations: 



(bias)p = 



(14a) 



where 

$Ma, P P) - const ■ exp (— (X> - P Ma ,dir(p)) 2 /24 p ) (15) 
so that 

(bias)p = r>Ma,dir(p) - Anv(p) (14b) 

This simple expression is then easy to integrate over all p in the 
sample, when one notes that the distribution of p is 



const • exp (— (p — p ) 2 / (2a' 2 0j>)) 



(16) 



where p' gives the Malmquist-shifted average value of p 
(shifted by a' times the shift of Do): 



p' Q = a'V + b' + (3 - a) ■ In 10 • a\ p 
Averaging (bias) p over all p from p\ to p u gives: 

(bias) avc = - (o Vp /(j v ) 2 ' avy/2/TrB(xi,x u ) 
where x\ and x u are 

Pl,u - Po 



V2a' 



(17) 



(18) 



(19) 



Note that the information on the space density (via a and Eq. 
16) is contained in x\ and x u , which show how far away from 
the Malmquist-shifted average p' the cutoffs are, in terms of the 
dispersion op = a'op of the cosmic (Gaussian) distribution 
function of p. The factor before B(x\, x u ) can also be written 
in terms of a p and a p. Then 

(bias)avc = -{al/itf + a 2 p)){a P /a')y/2/^B(x h x u ) (20) 

Such an analytic calculation is possible only in the given spe- 
cial cases of space density. Generally, the distribution of points 
in the p vs. loge? n diagram will differ from the Gaussian bi- 
variate case: at each p, the shift of (2?dir(p)) will be due to 
different Malmquist biases of the 1st kind. However, even in 
this case the inverse relation is not affected and the calculation 
and prediction of Sect. 3 is valid. 



Table 1. Bias in (log H) due to p cutoffs. The upper sign in the 
± symbols refers to p u while the lower sign refers to p\. 



Zu(,l) 


% 


-y/2/nB(-oo,Xu) 


bias* 






(-s/2/^B(x u oo)) 




±oo 


100 


0.00 


0.000 


±0.84 


80 


±0.22 


±0.015 


0.00 


50 


±0.80 


±0.056 


=F0.84 


20 


±1.68 


±0.118 



* for ax> = 0.28, cpd„ = 0.14. 



6. Application of the average bias formula 

We are interested in seeing how large an effect on log H cut- 
offs of different severities have. Table 1 shows the values of 
B(—oo,x u ) (B(x\,oo)) and the bias in (log if) for different 
values of upper (lower) cutoff ^ u (,l)i an d tne corresponding 
completeness percentage of the sample. 

In particular, if half of the galaxies are lost because of 
the upper cutoff (while there is no lower cutoff), formula (18) 
adopts a simple form, which gives a rough estimate of the bias 
for such distorted samples. 



(bias) 50% = -{a^Ja T ,) 2 a T )\f 7 2^B{~co^) 
= -0.8{<7v p /ctt>) 2 o-t> 



(21) 



With representative values ax> = 0.28 and a-p — 0.14 this 
gives the bias —0.056, or distances underestimated by about 
14 percent, in the average. Intuitively, one might have expected 
a more substantial influence of the p cutoff, knowing that the 
assumption of no p-selection is so fundamental for the inverse 
relation method. 



7. Problems of slope and calibration 

In the discussion thus far, we have assumed that the normal- 
ized distance has been accurately calculated, so that there are 
no observational errors in log D25 and the kinematical distance 
scale (dkin) is accurate as well. However, this cannot be the 
case, and as a result, field galaxies collected from a fixed true 
distance will show in the p vs. log d n diagram a different slope 
than without such errors. This slope will differ from the intrin- 
sic slope followed by calibrators. 

Two questions arise: first, in such a situation, what hap- 
pens to the slope-criterion requiring that (log H) for the unbi- 
ased plateau galaxies goes horizontally in the logiJ vs. dkin 
diagram? Secondly, how do we calculate the actual value of 
log Hq with calibrators known to have a different iTF slope? 
Her e we e xpand somewhat the compact treatment in Teeriko- 
rpi ( |l990t ). 

Consider two galaxy clusters in the p vs. log D25 diagram, 
both having (in their unbiased <i n -part), the same iTF slope a', 
which however differs from the calibrator slope. It is evident 
that the slope a' gives the correct relative distance of these clus- 
ters, when one shifts the regression lines one over the other (c.f. 
Fig. 1). Hence, as an answer to the first question, this is also 
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the slope which corresponds to the horizontal distribution of 
(log H); correct distance ratios reveal the linear Hubble law. 

Assume now that the calibrators follow the slope a' c > a'. 
One may describe in a picturesque manner how to do the 
comparison with calibrators, leading to an unbiased value of 
log Hq. As the shallower slope for the field galaxies is assumed 
to be due to errors in d n , one may think of "adding" random er- 
rors to calibrator linear log diameters V, until the slope falls to 
a' and then use the obtained relation for calibration ("shifting 
of regression lines"). Fortunately it is not necessary to actually 
"shuffle" 2?'s - this is equivalent to forcing the slope a' through 
the calibrators, at least in an ideal situation. 

So, let us assume first that the calibrator sample is volume- 
limited. Then adding random errors to T> means that 1(D) = T>$ 
remains the same, and one arrives at the relation 



Table 2. Hq as derived using a' = 0.5 when i?o(true) 
and the true calibrator slope is a' c . 



55 



(p) = a'{V) + b> 
On the other hand, the actual calibrators follow 

(j>)=a' c (V) + b' c 



(22) 



(23) 



Hence, after this trick, the new zero-point which must be ap- 
plied for the sample galaxies is: 



b' = b' c +(a' c -a'){V) 



(24) 



There is also a longer route to this result, using Eq. (7) giving 
the direct relation in terms of the inverse parameters and T>q . If 
one increases a-p so that a' c — > a' and requires that the direct 
relation does not change, one also arrives at Eq. (24). 

On the other hand, the formula for b' is exactly what is ob- 
tained if one simply forces the slope a' through the calibrators. 
From this follows that a correct calibration is achieved by de- 
riving the zero-point by a force-fit of the slope a' through the 
calibrators. 

8. The calibrator sample bias when one uses relevant field 
slope a' 

In the above discussion it is noteworthy that one does not have 
to know the intrinsic slope a' c for the calibrators. However, 
the calibrator sample must fulfill a special condition. It must 
be equivalent to a volume-limited sample, or reflect the cos- 
mic Gaussian distribution, as was already noted in Teerikorpi 
(1990). From Eq. (23), it is clear that two calibrator samples 
with different (D) yield different zero-points b' for the same a' 
of the field galaxies, while there can exist only one such correct 
b' , the one corresponding to (D) = Vq or (p) = a'Vo + b' = 
Po- 

If (2?) of the calibrator sample deviates from the cosmic 
T>q, there is a systematical error in 6', resulting in an average 
error in the derived (log H ) : 



A (log if) = (a'Ja'- 1)AV 
or in terms of the average p: 

A(logff) = {a'Ja'- l)Ap/a' 



(25) 



(26) 



1 

a c 


H {a = 0) 


H {a = 0.8) 


0.60 


67 


64 


0.65 


74 


68 


0.70 


82 


74 


0.75 


90 


79 



This source of error is very significant if the calibrator sam- 
ple has not been constructed with the intention of reaching the 
cosmic distribution of p, if for example the selection of galaxies 
has aimed at detection of Cepheids. In this way, one may have a 
volume-limited sample for each p (c.f. Theureau et al. 1997a), 
but this does not guarantee that the distribution of p reflects 
the cosmic distribution from which the field galaxy sample has 
been drawn. 

That this could be a genuine problem, consider the calibra- 
tor sample of KLUN with Cepheid distances. An indication of 
the calibrator sample bias is seen by comparing the distribution 
of calibrator p's with that of the whole sample. Because the lat- 
ter is diameter limited, it suffers from a Malmquist bias, hence 
its median should be displaced to a larg er valu e of p in compari- 
son with the calibrators. Ekholm et al. d!998| ) show that there is 
no such shift, which implies that calibrators' (p) is larger than 
the cosmic value. Such a situation results in an overestimated 
log H when one is compelled to use a' < a' c . 

Application of Eqs. (25-26) requires a knowledge of the 
average T>q or po of the cosmic distribution functions to be 
compared with those of the calibrator sample. It is clear that 
any calculation of T)q from the field sample requires, besides a 
suitable method, the value of Hq itself, hence an iterative ap- 
proach. Below we show that there is another much more con- 
venient route which does not require an explicit knowledge of 
the difference AT> (or the value of Hq). However, it assumes a 
radial space distribution law of galaxy number density. 

9. Another route: using calibrator slope a' c and making 
a' c /a' and Malmquist corrections 

In principle, the problem of the p-distribution of the calibrators 
does not affect the determination of the zero-point b' c when one 
forces the correct slope a' c on the calibrators. This makes one 
ask what happens if a' c and b' c are then used for the field sample. 
First of all one expects a dista nce-dependent bias, according to 
Eq.15 from Teerikorpi ( 1990| ), here adopted to diameters: 



A(logH)(r) = ({V)(r)-V )(a'Ja'-l) 



(27) 



Here (D) (r) is the average (log) diameter for the galaxies at 
true distance r. Now one may ask what is the average bias for 
the sample? The average bias clearly depends on the average 
of (D)(r) — Vq. This is actually the Malmquist bias of the 1st 
kind for Gaussian distribution of diameters with dispersion o~x>, 
and in the case of a space density distribution oc r~ a one gets: 

(A log If) = (3 -a) •ml0-o-f>(t4/a'- 1) (28) 



P. Teerikorpi et al.: Inverse Tully-Fisher relation 



7 



When a = the distribution is homogeneous. In order to use 
this as a correction, one must know both a' c and a! . Clearly, 
if the field and calibrator slopes are identical, then (log H) is 
unbiased, and this is so irrespective of the nature of the calibra- 
tor sample. Note also that this formula contains the width a-p 
of the total diameter function, not the smaller dispersion a-p ■ 
Hence the effect can be formidable even when the slopes a' c 
and a' do not differ very much. 

As a concrete example, consider the KLUN field galax- 
ies following the inverse slope a' as 0.5 which is in agree- 



ment with the simple disk model of Theureau et al. (1997a) 



This slope, when forced through the calibrators, gives Ho « 
75 kms" 1 Mpc" 1 . Above it was noted that the calibrator sam- 
ple cannot represent the cosmic distribution of log Vm, and this 
is the situation when a too high Ho is derived if the actual slope 
is steeper than the slope from the field galaxies. It is natural to 
ask, how large a slope a' c would in this manner explain the dif- 
ference between i/o(dir) w 55 and Ho(iav) w 75. Taking 
<Tj) as 0.25, we get predicted iJrj(inv) for different slopes a' c as 
given in Table 2. 

The last column corresponds to a decrease in the average 
space density, cor respo nding to density oc 7-~ - 8 as derived by 
Teerikorpi et al. ( 1998 ). It may be concluded that in order to 
explain the difference between Hq = 75 and 55 solely in this 
manner, it is required t hat ac tually a' c as 0.65 - 0.75, instead of 
ps 0.5. Ekholm et al. ( 1998[ ) show that such a steep calibrator 
slope really is consistent with the data. 



10. Concluding remarks 

In the present study we have searched for answers to a few ba- 
sic questions concerning the inverse TF relation as a distance 
indicator. First, how to see the signature of observational in- 
completeness in the TF parameter log Vm, analogously to the 
well-known increase in log H at increasing (normalized) dis- 
tance, in the case of the direct TF relation? Secondly, in which 
sense and how significant is the influence of the log Vm incom- 
pleteness alone on the derived value of log H (supposing one 
knows the correct inverse TF slope)? Thirdly, how to work with 
the situation when the inverse slope obeyed by the calibration 
sample differs from the relevant slope for the field sample (or 
for a cluster), and how important is in this regard the nature of 
the calibration sample? 

We have generalized the "fine-tuning" scheme discussed by 



ET97, and have pointed out that there is a useful log H vs. d n 



representation for the inverse relation, revealing both the rele- 
vant inverse slope and cutoffs in log Vm- Here d n is the normal- 
ized distance applicable to the inverse TF relation, actually the 
(log) linear diameter within an unknown constant. 

We have derived simple analytical formulae which give the 
influence of the logVM cutoffs on derived average distance 
(hence, on logH). The treatment is based on Gaussian distri- 
butions, and, in the case of the average bias, on the assumption 
that the space density is proportional to r~ a . The probable case 
that there is an upper observational cutoff in logVM, would 
cause an overestimated value of Hq. However, this effect is not 



at all dramatic - even if the upper cutoff excludes half of the 
galaxies, Ho would be increased typically by only 14 percent. It 
may be concluded that a cutoff in log Vm cannot alone explain 
derived large values of Ho- 

The problems of the relevant inverse TF slope and the na- 
ture of the calibrator sample lead to serious consequences on 
the use of the inverse TF relation as a distance indicator. It 
has been previously realized that the inverse slope a' applicable 
to the field sample is not necessarily the same as the slope a' c 
obeyed by the calibrators, which may result in a biased value 
of Ho (Teerikorpi 1990). In the present paper, such situations 



were systematically studied. 

If the two slopes are the same (in practice, when accuracies 
of calibrator distances are high and the photometric measure- 
ments for the field sample are not less accurate than those for 
the calibrators), then there is no bias, and it is also true that the 
nature of the calibrator sample (volume-limited, magnitude- 
limited, etc.) is not important ( ET97| ). 

If the sample inverse slope a' is shallower than the cali- 
brator slope (as normally expected), then matters come to de- 
pend critically on the nature of the calibrator sample. First, 
there is the possibility that the calibrator sample is a true, 
volume-limited representation of the Gaussian cosmic distribu- 
tion function of log Vm from which the field sample has been 
taken (the latter subject only to Malmquist bias). In this ideal 
case, it is permitted to use the field sample slope, with the zero- 
point obtained by a force-fit through the ca librators. This was 
the solution suggested by Teerikorpi ( 199C ), and in this case it 
is not necessary to know the value of the calibrator slope. 

However, generally the calibrator sample's (p) is shifted 
from the cosmic pq. If one knows the shift Ap (or AT>) and the 
slopes a' c and a', then one may calculate the resulting error in 
log-ff, when the field sample's a' is used. The systematic error 
is given by Eq. (25). 

Another, more convenient route, which bypasses the ex- 
plicit need for Ap, is to use the calibrators' slope a' c (and corre- 
sponding zero-point) and to calculate analytically the resulting 
systematic error in terms of the Malmquist bias of the 1st kind. 
This assumes that the average space number density around us 
is proportional to r~ Q . The systematic error is given by Eq. 
(27). 

We repeat that in the general case both a' c and a 1 (with cor- 
responding zero-points from force-fits through calibrators) re- 
sult in erroneous average distance estimates. Fortunately, we 
have now some quantitative control on this error. 
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